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Abstract 

The distance-resolved effective interaction between two star polymers 
in a good solvent is calculated by Molecular Dynamics computer simula- 
tions. The results are compared with a pair potential proposed recently 
by Likos et al. [Phys. Rev. Lett. 1998, 80, 4450] which is exponentially 
decaying for large distances and crosses over, at the corona diameter of the 
star, to an ultrasoft logarithmic repulsion for small distances. Excellent 
agreement is found in a broad range of star arm numbers. 
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Star polymers are hybrids between polymer-like entities and colloidal particles 
establishing an important link between these different systems, for recent reviews 
see ref |l],|^. The interpenetrability of two stars is governed by the number of 
arms (or functionality) /, i.e. the number of linear polymer chains attached to 
a central microscopic core. For / = 1,2 one recovers a system composed only 
of linear chains while in the limit f oo one gets sterically-stabilized spherical 
colloidal particles which behave like effective hard spheres.iB Recent researchi^i 
has mainly focussed on polymer conformations of a single star. The only relevant 
length scale of a single star is embodied in the spatial extension of the monomers 
around the core as given by the so-called corona diameter a. 

In order to predict macroscopic properties of a concentrated solution of many 
stars, one has, however, to proceed one step further: In any statistical theory, 
the effective interaction between the stars is a necessary input. This interaction, 
in general, comprises many-body terms. For concentrations which are not too 
high, i.e. smaller than or comparable to the overlap concentration p* = l/o"^, 
triplet and higher-order terms are small and the system is dominated by effec- 
tive pairwise interactions. Recently, based on scaling theory,!! an explicit ana- 



lytical expression for the effective pair potential V{r) was proposed in ref |T0 
This potential combines a logarithmic form of the interaction potential for core- 
core separations r smaller than a with an exponentially decaying interaction of 
Yukawa-form for distances r larger than a: 
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Note that the potential strength simply scales with the thermal energy ksT since 
the repulsion between the stars is of purely entropic origin having a good solvent 
in mind. Both, the potential in eq[^ and its associated force F{r) = —dV{r)/dr 
are continuous at r = cr, but F{r) has an artificial cusp at r = cr. The pref actor 
of the logarithm is fixed by scaling theory! while the exponential decay length 
la I ^ is the diameter of the outermost blob within the Daoud-Cotton model for 
one star polymer.! 

For an arm number of / = 18, this potential was tested against neutron 
scattering data and reasonable agreement was found.0 Further experimental 
support comes from shear moduli measurements in the crystalline phase of many- 
arm-micelles. Still, the scaling theory assumptions are strictly speaking only 
justified for core-core distances r much smaller than a, and the exponential decay 
length of the outermost blob size is an heuristic assumption. Hence the validity of 
the potential for arbitrary arm numbers can be questioned. In this paper, we test 
the pair potential against a microscopic model, resolving the monomers of the 
chains, by extensive Molecular Dynamics computer simulations. To be specific, 
we use a simulation model for star polymers developed by Grest et al.,i which was 
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applied in previous studies for single stars, and generalize it to a situation with 
two stars, which is the minimal set-up to extract information about the effective 
interaction between two stars. The distance-resolved interaction force F{r) is 
calculated for arm numbers / ranging from / = 5 to / = 50. Each arm contains 
N monomers where is varied from 50 to 200. As a result, we confirm the phe- 
nomenological interaction potential in eq |I|; our simulation results are in perfect 
quantitative agreement with the theoretical prediction. This important result en- 
ables a mapping of a star polymer solution onto a classical one-component fluid@ 
interacting via the effective ultra-soft pair potential of eq 0, provided the star 
concentration does not exceed p*. This picture was anticipated in recent work, 
calculating the anomalous structure factor of star polymer solutions^ and the 
unusual phase diagram including re-entrant melting^EJ and anisotropic crystal 
structures.0 So, our present work provides a theoretical justification of all these 
previous studies. 

Let us first describe the simulation model:! Each polymer arm consists of 
effective monomers or "beads" interacting via a purely repulsive Lennard- Jones- 
like potential Vo(r), where r is the separation of the beads. Vo(r) is obtained 
from the usual Lennard- Jones potential Vlj^t) by cutting Vi,j(r) at the position 
of the potential minimum = 2^/^(Tlj and by shifting it by the constant value 
VLj{rm) in order to obtain Vo(rm,) = 0: 




; r > 2i/Vij 



Here, e sets the energy scale and a^j the length scale of the beads. The pure 
repulsion implies that we are dealing with a good solvent. For neighbouring 
beads along the chains, the attractive FENE-potentialil Vch{r) is added to the 
interaction 
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This interaction diverges at r = i?o, which determines the maximum relative 
displacement of two neighbouring beads. Henceforth, we fix Rq to be 1.5aLj- 
Then the total potential Vo(r) + V^/i(r) between neighbouring monomers has a 
minimum at r ^ 0.97aLj. Furthermore, the central core particles of the two stars 
have a finite hard core radius Rc, and all monomers are interacting with the core 
particles via a modified repulsive interaction potential V^{r). The introduction 
of a small hard core of the central particles of the stars is necessary to accommo- 
date the large number of arms at small distances from the core.! Thus we take 
explicitly for the potential 

(4) 
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In addition, the innermost monomers of each arm are interacting with their core 
via an attractive potential which is given by 

We note that exactly this simulation model was already used by Grest et al. in 
their simulations of single star polymers.! In our simulations, the centers of the 
two stars are fixed at positions Ri and R2 with a given distance r = Ri — R2 . 
The total number of mobile monomers is which limits our studies to small 

/ and small A^. In all simulations, the system is held at fixed temperature 
T = 1.2e/ kB- Under these circumstances, the effective force Fi acting on the 
ith star center is given as a canonical average 

^. = (-v^:. (|EK,^(ki^-^.l) + EK^/.(l^i-^^l)j)' (6) 

where in the first sum the repulsive interactions of the core with all 2fN monomers 
in the system are considered, whereas the second sum only accounts for the at- 
tractive interactions with the / innermost monomers of the chains attached to the 
ith center. Obviously^due to symmetry, Fi = —F2. We use standard Molecular 
Dynamics simulationsH^ to equilibrate the monomers and perform the statistical 
average (...) over the monomers for the forces on the star centers. The timestep 



is typically 6t = 0.002r (with r = ^mcrlj/e being the associated time unit and 
m the monomer mass) and typically 120000 steps are used for equilibration and 
up to tjnax/St = 500000 steps were simulated to gather statistics. It was carefully 
checked by monitoring the internal energy that the system had equilibrated. A 
typical snapshot of two stars after equilibration is shown in Figure 0. As can be 
seen, the monomers of one star do not penetrate much into the central region of 
the other star. 

In order to check the code, we performed simulations of single stars changing 
the arm numbers between / = 5 and / = 50 and the monomer numbers from 
A^ = 50 to A^ = 200. The corresponding results for the radius of gyration 
-^G ~ jNJ2i=i{'^i ~ 'f^chiY (where tcm is the center of mass of the whole star) 
and the density profile of the monomers are in very good agreement with the 
results given in ref ^ and are well described by the scaling theory of Daoud and 
Cotton.El For a detailed list of the simulation parameters and the results for Rc^ 
obtained from these single star simulations, see Table |. 

It should be noted that the effective forces on the star centers are the gradient 
of the effective star-star potential. This effective potential, however, differs in 
general from the monomer averaged potential energy of the star centers. § We 
therefore had to calculate the averaged forces Fi {i = 1, 2) from our two star 
simulations to compare with the theoretical force as calculated from eq |I|. In 
doing this, two difficulties are arising: i) The corona diameter a, which is the 
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relevant length scale in the potential of eq 0, is not known a priori, ii) In contrast 
to the theory, there is a finite core size Rc in our simulation model. 

As regards the first difficulty, a is usually defined as the typical maximum 
range where a scaling behaviour of the monomer density around a single star 
center holds.Bi A statistical definition of a, however, is missing. On the other 
hand, the radius of gyration Rq has a clear definition as a canonical average, 
which can be calculated directly in simulations. We therefore use Rq, which was 
calculated in the single star simulations, as basic length scale for our simulation 
data and fit these data for F = Fi {i = 1, 2) to the theoretical prediction for 
F(r) using the least-square method and treating a as the single fit parameter. 
Afterwards we check how the optimal value for a scales with Rq as obtained from 
the single star simulations. The procedure is consistent if the ratio A = a /2Rg 
is independent of /. The second difficulty is resolved as follows: A logarithmic 
potential for r < a implies that the data should fall on a straight line crossing 
the origin if one plots the inverse force, 1/F, versus r inside the corona diameter. 
A typical plot is given in Figure ^. In fact, the data fall on a straight line. 
Extrapolating the data, however, one does not hit the origin. The divergence of 
the force occurs already at a finite distance 2Rd ~ 2Rc + (Jlj which clearly has to 
be attributed to the presence of the finite core in the simulations. We note that 
both, Rc and i?^, are microscopic length scales and of same order of magnitude 
(see Table D, thus being not relevant for the macroscopic length a in the scaling 
regime. We therefore normalize our distances by subtracting 2Rd, thus matching 
the divergence of the force properly. We emphasize that the slope of the straight 
line is in very good agreement with the theoretical prediction, see again Figure ^ 
This implies that the theoretical prefactor ^/^/^fc^T in eq |l] is confirmed by the 
computer simulations. In Figure ^, we show the effective force versus distance for 
five different arm numbers / and two monomer numbers N . The agreement with 
the theory is convincing for all / and A^. The consistency of our fitting procedure 
of the corona diameter a is documented in Table |, where the ratio A = a I2Rq is 
given for different / and N . We find A ~ 0.65 independent of /. This value also 
coincides with the value used in ref [1^ to fit experimental data for / = 18. We 
further note that A is independent of A^, consistent with scaling theory. Finally, 
we prove the exponential decay of the force for distances larger than a by plotting 
the logarithm of the force versus distance in Figure ^ for one typical example. 
One clearly sees the crossover of the inner-core data to a straight line outside 
the core. The slope is consistent with the theoretical one as determined by the 
outermost blob size. 

In conclusion, we have verified the ultra-soft pair interaction for star polymers 
by direct molecular simulations. It is straightforward to generalize the method 
to two stars confined in a periodically-repeated cubic cell in order to estimate 
the shrinking of the corona diameter due to a finite star density. Also, similar 
to charged colloids,0 triplets of stars should be considered to investigate the 
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importance of triplet interactions. Our future work lies along these directions. 

Wc thank M. Schmidt and C. N. Likes for helpful remarks. 
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Tables 



/ 


N 


6t/T 


tmax/ ^t 


Rg/ctlj 


{Rc/<Jlj) + 1/2 


Rd/CTLJ 


A 


5 


100 


0.004 


500000 


13.53 


0.65 


1.39 


0.61 


10 


50 


0.003 


400000 


10.37 


1.1 


1.21 


0.66 


10 


100 


0.003 


400000 


16.18 


1.1 


0.89 


0.64 


10 


150 


0.003 


400000 


19.71 


1.1 


1.31 


0.60 


10 


200 


0.003 


400000 


24.52 


1.1 


1.42 


0.67 


18 


50 


0.002 


350000 


11.19 


1.25 


1.38 


0.68 


18 


100 


0.002 


350000 


17.10 


1.25 


1.64 


0.65 


30 


50 


0.002 


350000 


12.22 


1.6 


1.89 


0.66 


50 


50 


0.002 


350000 


13.35 


1.8 


2.40 


0.69 



Table I: List of the simulation parameters and the corresponding results for Rq 
and A = a/^Ra- 



Figure captions 

Figure 1: Typical configuration for two stars with / = 10 and = 50. The 
distance between the central core particles, which are shown as big black spheres, 
is r = 5.2(Jlj. The gray and light gray monomers are belonging to the first and 
second star respectively. 



Figure 2: Reduced inverse force ksT / {FRq) between the centers of two star 
polymers (for / = 10 and N = 50) versus reduced distance r / Rq. The error bars 
were obtained by averaging over the results of 8 independent simulations. 



Figure 3: Simulation results (symbols) and theoretical results (lines) for the 

reduced effective force FRc/ksT versus reduced distance (r — 2Rd)/RQ. a) for 
/ = 5, 10, 18 and N = 100, b) for / = 18, 30, 50 and N = 50. 



Figure 4: Logarithm of the reduced force hi. (FRc/kBT) versus reduced dis- 
tance (r — 2Ra)/RG for / = 10 and N — 50. The error bars were obtained by 
averaging over the results of 8 independent simulations. 
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Jusufi et al., fig.2 
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Jusufi et al., fig.Sa 




Jusufi et al., fig.Sb 




Jusufi et al., fig.4 
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